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Two  important  observations  when  the  onset  of  spacecraft  surface  charging  occurs  arc  (1)  the 
electron  flux  measured  in  the  high-energy  (above  several  kcV)  channels  onboard  rises  and  (2)  the 
ambient  electron  temperature  rises  above  a  critical  value.  We  show  by  means  of  an  analytical  model 
that  the  two  behaviors  arc  consistent  with  each  other.  [DOI:  10.1063/1.3125517] 


I.  INTRODUCTION 

It  has  been  observed  repeatedly  on  the  Spacecraft  Charg¬ 
ing  at  High  Altitudes  (SCATHA)  satellite  that  when  space¬ 
craft  surface  charging  occurred,  the  electron  fluxes  measured 
in  the  channels  of  above  several  kcV  rose  while  those  mea¬ 
sured  in  the  low  energy  channels  below  several  kcV  fell.  In 
cases  of  severe  charging,  the  rise  in  flux  obtained  in  the 
channels  of  about  20-30  keV  was  remarkable  and  had  re¬ 
ceived  much  attention.  These  observations  have  led  to  some 
widespread  belief  that  a  rise  in  the  flux  of  the  20-30  keV 
channel  is  one  of  the  best  symptoms  to  signal  the  occurrence 
of  spacecraft  charging.  On  the  other  hand,  one  can  show 
analytically  by  means  of  the  Maxwellian  model  that 
spacecraft  charging  occurs  when  the  ambient  electron  tem¬ 
perature  rises  above  a  critical  temperature  and  that  the  charg¬ 
ing  level  increases  as  the  temperature  rises  further.  Observa¬ 
tions  on  the  Los  Alamos  National  Laboratory  (LANL) 
geosynchronous  satellites  have  provided  abundant  evidences 
supporting  the  temperature  dependence  idea  of  spacecraft 
charging.  Each  of  the  two  behaviors  mentioned  above  is 
important,  but,  like  the  Indian  fable  of  four  blind  persons 
touching  an  elephant,  the  two  behaviors  seem  to  be  unrelated 
aspects  in  spacecraft  charging.  Some  even  thought  that  the 
two  ideas  arc  mutually  exclusive.  The  purpose  of  this  note  is 
to  show  in  a  theoretical  framework  that  the  two  behaviors  arc 
mutually  consistent.  In  the  next  section,  we  begin  with  the 
definition  of  flux,  develop  the  analytical  formulation,  and 
present  numerical  solutions  and  their  physical  interpreta¬ 
tions. 


II.  MAXWELLIAN  FLUX 

The  flux  J  is  defined  in  terms  of  the  Maxwellian  velocity 
distribution  F(v)  of  the  electrons, 

J  =  j  d'vvFiv).  (1) 

In  Eq.  (1),  the  velocity  distribution  F(v)  is  given  by 

a)E!ectronie  mail:  afrI.rvb.pa@hanseom.af.mil. 


F(v)  =  n 


m 

2vkT 


(2) 


The  velocity  v  can  be  expressed  in  terms  of  energy  E  with 
the  substitution  £=(1 I2)mv2,  After  the  substitution,  Eq.  (2) 
is  converted  to  the  Maxwellian-Bolt/mann  distribution  f(E) 
of  the  ambient  electrons. 


f{E)  =  n 


iSr)  cxp(-s)- 


(3) 


where  T  is  the  electron  temperature,  m  the  electron  mass,  and 
E  the  electron  energy.  Figure  1  shows  a  plot  of  lo g,/(£) 
versus  E.  The  plot  shows  a  straight  line  and  its  slope  equals 
-1/AT.  The  y  intercept  at  £=0  is  given  by 


/«>)  =  » 


2irkT 


(4) 


which  is  a  function  of  the  density  ?t  and  the  Maxwellian 
temperature  T  of  the  ambient  electrons.  As  the  temperature  T 
changes,  the  location  of  the  y  intercept  changes.  Note  that 
one  can  use  log10/(£)  and  the  plot  would  be  a  straight  line. 


FIG.  I.  (Color  online)  Log  of  the  Maxwellian  energy  distribution  function 
plotted  against  energy  E.  The  slope  of  the  straight  line  depends  on  the 
temperature  T.  The  y  intercept  is  a  function  of  the  density  n  and  the  tem¬ 
perature  T.  A  plot  of  log  K{E)  vs  E  is  also  plotted.  It  features  a  maximum  at 
E-T .  (For  simplicity,  both  n  and  T equal  1  in  the  input.  To  show  the  eorreet 
slope,  one  needs  to  use  the  natural  log  seale  ) 
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but  there  is  a  multiplicative  constant  because  lo geJ'(E) 
=  logf  10Xlog1()/(£). 

For  normal  incidence,  the  flux  J  of  Eq.  (1)  can  be  sim¬ 
plified.  Neglecting  the  multiplicative  constant  of  it  and  /w, 
which  do  not  vary  during  spacecraft  charging,  we  write  the 
normal  flux  J  [Eq.  (3)]  as  a  function  of  /(E), 

J=  \  dEEf(E).  (5) 

Jo 

The  flux  J ,  Eq.  (5),  can  be  written  as 


dE 


dJ(E) 
dE  ' 


(6) 


where  the  kernel  ( dJIdE )  of  the  integral,  Eq.  (5),  is  called 
the  differential  flux  K(E ), 


K(E)  = 


dJ(E) 

dE 


=  Ef(E ). 


(7) 


A  plot  of  log  K(E)  is  shown  in  Fig.  1. 


III.  MEASUREMENT  TECHNIQUE 

In  practice,  the  instruments  for  flux  measurement  on  sat¬ 
ellites  have  limitations  Instead  of  obtaining  a  continuous 
distribution  of  flux  J{E)  as  a  function  of  energy,  the  measure¬ 
ment  obtains  the  flux  in  discrete  energy  channels,  each  chan¬ 
nel  having  its  finite  energy  range  A E.  That  is,  the  measure¬ 
ment  gives  A7(E)/AE,  where  A  represents  a  finite  increment. 
We  can  rewrite  Eq  (5)  as  follows: 

■/  =  2  ( =  2  K£E)lEh  (8) 

where  [A7(E)/AE]r  denoted  by  E,(E),  is  the  differential  flux 
in  the  ith  channel, 


An  example  of  differential  flux  measurements  is  shown  in 
Fig.  2.  In  this  example,  the  ambient  flux  is  measured  in  dis¬ 
crete  energy  channels.  The  ith  channel  measures  the  differ¬ 
ential  flux  K,(E)  at  Et  ±AEj. 


IV.  DILEMMA 

Each  of  the  two  techniques  for  signaling  the  occurrence 
of  spacecraft  charging  has  received  credible  observational 
evidence.  In  the  high  energy  channel  technique,  when  a  high 
energy  channel,  such  as  20-30  keV,  rises  relative  to  the 
lower  energy  channels  (keV  or  less),  spacecraft  charging  oc¬ 
curs  (Fig.  2).  In  the  critical  temperature  technique,  when  the 
ambient  electron  temperature  rises  above  a  critical  value, 
charging  occurs  (Fig.  3).  We  will  not  attempt  to  simulate 
case-specific  results  because  more  parameters  are  often 
needed.  Wc  point  out  that  both  techniques  have  received  con¬ 
crete  evidence  obtained  in  space  and  yet  they  seem  to  be 
unrelated  to  each  other.  There  is  need  to  bridge  an  under¬ 
standing  concerning  both  techniques. 


While  there  is  a  Maxwellian  model'  for  explaining  the 
physics  in  the  “critical  temperature"  technique,  the  high  en¬ 
ergy  channel  technique  is  based  solely  on  experimentalist's 
observations.  We  will  attempt  to  prove  the  mutual  consis¬ 
tency  of  the  two  techniques  in  the  following  sections. 

Before  we  discuss  flux,  let  us  examine  the  relation  of  a 
Maxwellian  distribution  /(E),  Eq.  (1),  with  temperature. 
Since  the  slope,  d[\ogef(E)]/dE,  equals  -1/AT,  rising  tem¬ 
perature  would  correspond  to  higher  population  at  high  en¬ 
ergies  and  lower  population  at  low  energies  [Fig.  4(a)].  For 
two  temperatures  Tj  and  T2(>T j),  the  transition  energy  Et ,  at 
which  the  distribution  /(£,)  is  unchanged,  is  given  by 


(3/2)log  /7VT,) 
(1/AT|)  -  (1/AT2) 


(10) 


If  the  distribution  deviates  from  Maxwellian  and  one  can 
still  fit  an  average  slope  to  the  curve  of  /(E),  the  “critical 
temperature"  technique  is  still  approximately  valid.  If  the 
distribution  deviates  in  an  unmanageable  manner,  tempera¬ 
ture  is  undefined  [Fig.  4(b)]  and  a  manageable  way  to  char¬ 
acterize  the  electron  distribution  may  be  impossible. 


V.  EFFECT  OF  RISING  TEMPERATURE  ON  THE  NET 
INCOMING  FLUX 


First,  to  examine  the  effect  of  rising  temperature  T  on 
the  net  incoming  flux  K(E)  in  a  channel,  wc  differentiate  Eq. 
(7)  by  T, 


rlK(E.T)  /  3  j  1 

- =  E -  -kT\——^K(E). 

<TT  \  2  /(kT)2 


(11) 


This  simple  result,  Eq.  (11),  reveals  an  interesting  behavior. 
The  temperature  gradient  of  differential  flux,  E(E),  is  nega¬ 
tive  for  E  below  (3/2)AEbut  positive  for  E  above  (3/2 )AT 
As  the  temperature  rises,  the  differential  flux  K  increases  for 
T  below  2E/3  but  decreases  for  T  above  2E/3.  At  kT 
=  2E/3,  the  differential  flux  is  maximum.  Figure  5  illustrates 
the  effect. 

Next,  to  seek  out  the  most  sensitive  channel  for  a  given 
temperature,  we  differentiate  K(E)  by  E, 


dK(E,T) 

rlE 


(12) 


Equation  (12)  reveals  that  the  differential  flux  increases  with 
the  energy  E  as  long  as  E  is  below  kT.  For  E  above  AT,  the 
differential  flux  decreases  as  the  energy  increases.  The  maxi¬ 
mum  of  the  differential  flux  K  is  at  E=kT.  Figure  7  illustrates 
the  effect. 

With  spacecraft  charging  to  a  potential  (/>  (negative  volt), 
the  incoming  electron  flux  is  reduced  by  the  Boltzmann  fac¬ 
tor,  exp(-qe</>/  AT), 

K{E.T)  =  £/(£)exp|-^p).  (13) 

Accordingly,  the  location  of  the  zero  gradient  of  Eq.  (11)  is 
shifted  by  qe</>,  where  qe  is  the  charge  of  electron.  That  is,  the 
location  of  Max[r?E/r?T]  is  given  by 
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FIG.  2.  Fluxes  measured  in  channels 
of  different  energies  before  and  after  a 
significant  charging  event  Charging 
occurs  at  6:30  to  about  7:()OUT.  The 
temperature  of  the  main  component 
Maxwellian  rises  from  9.2  keV 
(6:30UT)  to  23  keV  (6  50UT)  (from 
Ref.  1). 


UT 


3 

E  +  qe(j>--kT  =  0.  (14) 

The  location  of  Maxf dKfdE]  given  by  Eq.  (12)  is  unchanged 
by  the  potential  d>. 

VI.  RESULTS 

We  have  obtained  numerical  results  in  the  theoretical 
framework  of  Maxwellian  plasmas.  Figure  5  shows  the  cal¬ 
culated  differential  flux  K{E,T)  as  a  function  of  temperature 
T  for  various  energies  corresponding  to  typical  channel  en¬ 
ergy  values.  The  spacecraft  potential  is  assumed  to  be  zero. 
The  lowest  energy  channel  (0.6  keV)  shows  the  highest  flux 
and  the  higher  energy  channel  shows  lower  fluxes,  as  one 
expects  from  the  general  nature  of  Maxwellian  or  any  equi¬ 
librium  distributions.  At  higher  temperatures,  the  fluxes  at 


higher  energy  (£'=9.2,  23.8  keV)  channels  rise  whereas 
those  at  lower  energies  fall.  For  a  given  energy  £\  the  maxi¬ 
mum  of  the  flux  curve  is  located  at  T=(2E73),  in  agreement 
with  Eq.  (11). 

Figure  6  shows  the  same  curves  but  with  the  spacecraft 
potential  equal  to  -4  kV.  Physically,  the  (luxes  of  the  low 
energy  channels  (£=0.6,  1.6  keV)  arc  suppressed  clearly  be¬ 
cause  the  low  energy  ambient  electrons  can  hardly  reach  the 
spacecraft  that  is  charged.  Mathematically,  the  flux  reduction 
is  by  the  multiplicative  Boltzmann  factor.  The  maxima  of  the 
curves  arc  no  longer  at  T=(2E!3)  but  at  T=2{E+i{e(f>)/3,  in 
agreement  with  Eq.  (14).  A  prominent  feature  in  this  figure  is 
the  relative  enhancement  of  the  fluxes  in  the  higher  energy 
channels 

Figure  7  shows  the  (luxes  as  a  function  of  energy  for 
various  given  temperatures.  The  spacecraft  is  considered  to 
be  uncharged.  As  expected  from  Fig.  2,  the  higher  energy 
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Jun  6  Jun  8  Jun  10  Jun  12 


FIG.  3.  (Color  online)  Charging  on  (he 
famous  Bastille  Day  event.  2000.  The 
spacecraft  potential  (negative  kV) 
rises  whenever  the  electron  tempera¬ 
ture  exceeds  a  critical  value  (from  Ref 

7). 


electron  population  increases  with  the  temperature  and  so  do 
the  high  energy  fluxes,  whereas  the  low  energy  fluxes  do 
oppositely.  The  maximum  of  each  curve  is  located  at  E=T, 
in  agreement  with  Eq.  (12).  Note  also  that  the  area  under  the 
curve  K(E)  is  given  by  the  following  integral: 


FIG.  4.  (Color  online)  (a)  L  og  of  the  Maxwellian  energy  distribution  func¬ 
tion  plotted  against  energy  E  for  two  different  temperatures  As  temperature 
increases,  both  the  intercept  and  the  slope  decrease.  The  function  increases 
for  E  exceeding  3*772.  (b)  Deviation  from  Maxwellian.  If  the  deviation  is 
gentle,  an  average  slope  can  still  be  defined.  The  idea  of  "higher  temperature 
increases  the  high  energy  population  ’  is  still  valid.  For  distributions  with 
unmanageable  deviations,  temperature  is  undefined 


(15) 


The  square  bracket  in  Eq.  (15)  is  unity.  Therefore,  the  inte¬ 
gral  of  Eq.  (15)  is  proportional  to  T12.  That  is,  the  flux  J 
increases  as  the  electron  temperature  T  increases,  as  ex¬ 
pected. 

Figure  8  shows  the  same  curves  as  in  Fig.  7  but  with  the 
spacecraft  potential  at  -4  kV.  Again,  the  suppression  of  the 
low  energy  fluxes  is  because  of  the  Boltzmann  factor.  The 
higher  energy  (£=9.2,23.8  keV)  fluxes  at  20-30  keV,  for 
example,  arc  enhanced  relatively. 


TEMPERATURE  (keV) 


FIG.  5.  (Color  online)  Flux  in  four  energy  channels  as  a  function  of  electron 
temperature  T.  The  spacecraft  is  uncharged.  As  the  energy  E  rises,  the  flux 
maximum  at  7=2£73  shifts  to  the  right  hand  side. 
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FIG.  6.  (Color  online)  Flux  in  four  energy  channels  as  a  function  of  elcetron 
temperature  T.  The  spacecraft  is  charged  lo  -4  kV  As  the  energy  E  rises, 
the  flux  maximum  shifts  to  higher  value  of  T. 

VII.  CONCLUSION 

From  this  study,  we  conclude  that  the  ideas  of  rising 
temperature  and  relative  enhancement  of  high  energy  fluxes 
arc  not  mutually  exclusive.  Even  without  spacecraft  charg¬ 
ing,  the  high  energy  channel  fluxes  increase  with  the  ambient 
electron  temperature.  Furthermore,  if  spacecraft  charging  oc¬ 
curs,  the  low  energy  channel  fluxes  are  suppressed  relative  to 
the  high  energy  channels  because  of  Coulomb  repulsion  ren¬ 
dering  the  enhancement  of  the  high  energy  channel  fluxes 


ENERGY  E  (keV) 


FIG  1,  (Color  online)  Differential  flux  as  a  function  of  energy  at  various 
temperatures.  The  spacecraft  is  uncharged  (0=0).  Note  that  the  maximum 
of  each  curv  e  is  at  E=kT. 
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FIG.  8.  (Color  online)  Same  caption  as  the  previous  figure,  hut  the  space¬ 
craft  potential  is  -4  kV.  The  amplitude  Is  reduced  by  exp kT).  Note 
lhat  the  maximum  of  each  curve  is  at  E-kT  and  the  flux  of  T=  0.6  keV  Is 
reduced  to  praelieally  zero. 

more  prominent.  Whereas  the  method  of  critical 
temperature  gives  precise  signaling  of  the  onset  of  space¬ 
craft  charging  as  verified  8  by  the  Los  Alamos  National 
Laboratory  (LANL)  geosynchronous  satellite  data,  the 
method  of  flux  enhancement  at  high  energies1  gives  promi¬ 
nent  signals  when  charging  occurs  but  without  the  precision 
of  the  onset  of  charging.  It  would  be  misleading  to  interpret 
from  the  rising  flux  in  the  20-30  keV  channels  during  space¬ 
craft  charging  that  only  the  20-30  keV  electrons  arc  respon¬ 
sible  for  charging.  It  is  advisable  to  “see  the  full  picture”  by 
considering  the  entire  electron  energy  distribution  for  mod¬ 
eling  and  forecasting  spacecraft  charging. 
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